Dephasing in the semiclassical limit is system— dependent 
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We investigate dephasing in open quantum chaotic systems in the limit of large system size to 
Fermi wavelength ratio, L/Af 3> 1. We semiclassically calculate the weak localization correction 
g wl to the conductance for a quantum dot coupled to (i) an external closed dot and (ii) a dephasing 
voltage probe. In addition to the universal algebraic suppression g wl oc (I + td/t^) -1 with the dwell 
time td through the cavity and the dephasing rate t7/ , we find an exponential suppression of weak 
localization by a factor oc exp[— r/r^], with a system-dependent f. In the dephasing probe model, 
f coincides with the Ehrenfest time, f oc lnfL/Ap], for both perfectly and partially transparent 
dot-lead couplings. In contrast, when dephasing occurs due to the coupling to an external dot, 
f oc ln[L/£] depends on the correlation length £ of the coupling potential instead of Af- 
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Introduction. Electronic transport in mcsoscopic 
systems exhibits a range of quantum coherent effects such 
as weak localization, universal conductance fluctuations 
and Aharonov-Bohm effects 0, 0] ■ Being intermediate in 
size between micro- and macroscopic systems, these sys- 
tems are ideal playgrounds to investigate the quantum- 
to-classical transition from a microscopic coherent world, 
where quantum interference effects prevail, to a macro- 
scopic classical world Indeed, the disappearance of 
quantum coherence in mesoscopic systems as dephasing 
processes set in has been the subject of intensive theo- 
retical [1, @, @, 0] and experimental @, 0, [13] studies. 
When the temperature is sufficiently low, it is accepted 
that the dominant processes of dephasing are electronic 
interactions. In disordered systems, dephasing due to 
electron-electron interactions is known to be well mod- 
eled by a classical noise potential Q , which gives an al- 
gebraic suppression of the weak localization correction to 
conductance through a diffusive quantum dot, 



.,wl 



= ffo wl /(l + r D /T ). 



(1) 



Here, g™ 1 is the weak localization correction (in units of 
2e 2 /h) without dephasing, the dephasing time is given 
by the noise power, and td is the electronic dwell time 
in the dot. Eq. ([1]) is insensitive to most noise-spectrum 
details, and holds for other noise sources such as electron- 
phonon interactions or external microwave fields. 

Other, mostly phenomcnological models of dephasing 
have been proposed to study dephasing in ballistic sys- 
tems @, H, 0] , the most popular of which, perhaps, being 
the dephasing lead model 0, @| . A cavity is connected 
to two external, L (left) and R (right) leads of widths 
Wl,Wr. A third lead of width W3 is connected to the 
system via a tunnel-barrier of transparency p. A voltage 
is applied to the third lead to ensure that no current flows 
through it on average. A random matrix theory (RMT) 
treatment of the dephasing lead model leads to Eq. (Q]) 
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Figure 1: a) (left panel) Schematic of the system-environment 
model. The system is an open quantum dot that is coupled 
to an environment in the shape of a second, closed quantum 
dot. b) (right panel) Schematic of the dephasing lead model. 



with td = t L/(Wl + Wr) and ts = t L/ (/0W3), in term 
of the dot's time of flight tq [fj], |Tl|. Thus it is commonly 
assumed that dephasing is system- independent. The de- 
phasing lead model is often used phenomenologically in 
contexts where the source of dephasing is unknown. 

Our purpose in this article is to revisit dephasing in 
open chaotic ballistic systems with a focus on whether 
dephasing remains system-independent in the semiclas- 
sical limit of large ratio L/Xp of the system size to 
Fermi wavelength. This regime sees the emergence of 
a finite Ehrenfest time scale, r^ 1 = A _1 ln[L/Ap] (A is 
the Lyapunov exponent), in which case dephasing can 
lead to an exponential suppression of weak localization, 
oc exp[— t^/t^]/ (1 + td/t^) [T^ . Subsequent numerical 
investigations on the dephasing lead model support this 
prediction Here we analytically investigate two dif- 
ferent models of dephasing, and show that the suppres- 
sion of weak localization is strongly system-dependent. 
First, we construct a new formalism that incorporates 
the coupling to external degrees of freedom into the scat- 
tering approach to transport. This approach is illustrated 
by a semiclassical calculation of weak localization in the 
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case of an environment modeled by a capacitively cou- 
pled, closed quantum dot. We restrict ourselves to the 
regime of pure dephasing, where the environment does 
not alter the classical dynamics of the system. Second, 
we provide the first semiclassical treatment of transport 
in the dephasing lead model. We show that in both cases, 
the weak localization correction to conductance is 



wl 



So 1 exp[-f/r^] /(l +t d /t ), 



(2) 



where is the finite- correction in absence of de- 
phasing. The time scale f is system-dependent. For the 
dephasing lead model, f = + (1 — p)r^ p in terms of the 
transparency p of the contacts to the leads, and the open 
system Ehrcnfest time t^ p = A -1 ln[VK 2 /ApL]. This an- 
alytic result fits the numerics of Ref. [IH , and (up to log- 
arithmic corrections) is in agreement with Ref. [121 ]. Yet 
for the system-environment model, f = A -1 ln[(L/£) 2 ] 
depends on the correlation length £ of the inter-dot cou- 
pling potential. We thus conclude that dephasing in the 
semiclassical limit is system-dependent. 

Transport theory for a system-environment 

model. In the standard theory of decoherence, one 
starts with the total density matrix r/tot including both 
system and environment degrees of freedom The 
time-evolution of ?7tot is unitary. The observed proper- 
ties of the system alone are given by the reduced density 
matrix ?7 sys , obtained by tracing ?7tot over the environ- 
ment degrees of freedom. This is probability conserving, 
Tr Vsys = 1, but renders the time-evolution of 7y sys non- 
unitary. The decoherence time is inferred from the decay 
rate of its off-diagonal matrix elements Q . We generalize 
this approach to the scattering theory of transport. 

To this end, we consider two coupled chaotic cavities 
as sketched in Fig. la. Few-electron double-dot systems 
similar to the one considered here have recently been the 
focus of intense experimental efforts [14J. One of them 



(the system) is an open quantum dot connected to two 
external leads. The other one (the environment) is a 
closed quantum dot, which we model using RMT. The 
two dots arc chaotic and capacitively coupled. In par- 
ticular, they do not exchange particles. We require that 
Af <C Wl,r <C L, so that the number of transport chan- 
nels satisfies 1 <C A^l,r <C L/Xf and the chaotic dynam- 
ics inside the dot has enough time to develop, Atd 3> 1. 
Electrons in the leads do not interact with the second dot. 
Inside each cavity the dynamics is generated by chaotic 
Hamiltonians H sys and H env . We only specify that the 
capacitative coupling potential, U, is smooth, and has 
magnitude U and correlation length £. 

The environment coupling can be straightforwardly in- 
cluded in the scattering approach, by writing the scatter- 
ing matrix, S, as an integral over time-evolution opera- 
tors. We then use a bipartite semiclassical propagator to 
write the matrix elements of § for given initial and final 
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environment positions, (qo,q), as 
§mn(qo,q) = (2tt) _1 / dt dy / dy (m|y) (y 

JO Jh JR 

x ^(C 7 C r )iexp[i{S , T + 5r + S 7 ,r}].(3) 

This is a double sum over classical paths, labeled 7 for 
the system and T for the environment. For pure dephas- 
ing, the classical path 7 (r) connecting y (qo) to y (q) 
in the time t is solely determined by H sys (H env ). The 
prefactor C 7 Cr is the inverse determinant of the stabil- 
ity matrix, and the exponent contains the non-interacting 
action integrals, S 7 , Sr, accumulated along 7 and T, and 
the interaction term, S 7 ,r = /QdT l lt[y.y(T),qr(T)]. 

Since we assume that particles in the leads do not inter- 
act with the second cavity, we can write the initial total 
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with Tysyi = \ri)(n\, 



density matrix as 
n = 1, 2, ...Nt,. We take r] env as a random matrix, though 
our approach is not restricted to that particular choice. 
We define the conductance matrix as the following trace 
over the environment degrees of freedom, 



M _ 

ztmn 



= (m|Tr cnv [S? ? ( ' l) S f ] \m). 



(4) 



The conductance 



is then given by g = Xi ra ,,J™- 
This construction is current conserving, however the 
environment-coupling generates decoherence and the 
suppression of coherent contributions to transport. To 
see this we now calculate the conductance to leading or- 
der in the weak localization correction. 

We insert Eq. ([3]) into Eq. ([4|, perform the sum 
over channel indices with the semiclassical approxima- 

tion En L (y ol n )("lyp) ~ s (y'o - yo) 0, and use the 

RMT result (q n l?wlqn) « O^nv^CQo _ Qo), where fi onv 

16[. The conductance then 



(qohcnvK) ~ 

is the environment volume 
reads 

g = (4^ 2 O cnv ) _1 f dtdt' f dq dq /dy f dy 

JO JSJ cnv Jh JR 

x J2 (C 7 C r C y CrO* e^^+^+^l (5) 

7T;7',r' 

This is a quadruple sum over classical paths of 
the system (7 and 7', going from yo to y) and 
the environment (r and V , going from qo to q), 
with action phases $ sys = S 7 (y ,y;t) - S 7 > (y , y; *'), 
$env = 5 r (qo,q;^) - 5r'(qo,q;i') and $ u = 
§ 7 ,r(yo,y;qo,q;*) - § 7 'T'(yo,y;qo,q;^')- We are in- 
terested in the conductance averaged over energy vari- 
ations, and hence look for contributions to Eq. ([5]) with 
stationary <& sys ,$ env . The first such contributions are 
the diagonal ones with 7 = 7' and r = V, for which 
$>u = 0. They are U-independent and give the clas- 
sical, Drude conductance, g° = Ni,Nr/(Nl + Nr). 
The leading order correction to this comes from weak- 
localization paths 7 and 7' [H, EH EJ El (see Fig. [2J), 



3 




Chaotic system (environment not shown) 



Figure 2: (color online) A semiclassical contribution to weak 
localization for the system-environment model. The paths are 
paired everywhere except at the encounter. There one crosses 
itself at angle e, while the other does not (going the opposite 
way around the loop). Here we show £ > eL, so the dephasing 
(dotted path segment) starts in the loop (T(_ > 0). 



with r' = r for environment paths [19j. In the absence 
of dephasing these contributions accumulate a phase dif- 
ference <5$ sys . Semiclassically these contributio ns g ive 

.9o w! = 



exp[-r^/r D ] N L N R /(N L + N R ) 2 0, 03 [T 



In the presence of an environment, each weak local- 
ization pair of paths accumulates an additional action 
phase difference <53>u, which is averaged over. Dephas- 
ing occurs mostly in the loop, when the paths are more 
than the correlation length £ apart (see Fig. [5]). Thus 
we can define 7^ = A -1 ln[(£/eL) 2 ] as twice the time 
between the encounter and the start of dephasing. If 
£ < eL, dephasing starts before the paths reach the 
encounter, T$ < [20]. Using the central limit theo- 
rem and assuming a fast decaying interaction correlator, 
(U[y 7 (r), q r (T)]U[y 7 (r'), qr(r')]) r oc t" 1 exp[-{y 7 (r) - 
y 7 (r' )}/£], the average phase difference due to 11 reads 



exp[-(t 2 -ii-T c )/7*], (6) 



where t\ fa) gives the start (end) of the loop. The deriva- 
tion then proceeds as for U = [ID], except that during 
the time fa — tx — T%), the dwell time is effectively di- 
vided by [1 + td/t^], so the fa — fi)-integral generates 
an extra prefactor of [1 + td/t^] -1 cxp[— r^/r^] where 
T£ = A -1 ln[(L/£) 2 ]. Thus the weak localization correc- 
tion is as in Eq. (|2|) with f given by t£. In contrast to 
Ref. [IH, the exponent depends on £ not Ap. 

A calculation of coherent-backscattering with dephas- 
ing to be presented elsewhere enables us to show that 
our approach is probability- and thus current-conserving. 
We also point out that (for ~ tb) one can ignore the 
modifications of the classical paths due to the coupling 
to the environment, as long as £ 3> [AfL/Atd] 1 / 2 . Thus 
our method is applicable for £ smaller (as well as larger) 
than the encounter size [ApL] 1 / 2 , but not for £ ~ Ap. 

Dephasing lead model. We next add a third lead 
to an otherwise closed dot (as in Fig. [TJa) , and tune the 
potential on this lead such that the net current through 



it is zero. Thus every electron that leaves through lead 
3 is replaced by one with an unrelated phase, leading to 
a loss of phase information without loss of particles. In 
this situation the conductance from L to R is given by 
9 = Tlr + TL 3 T R3 (T L3 + r R3 ) _1 @], where T nm is the 
conductance from lead m to lead n when we do not tune 
the potential on lead 3 to ensure zero current. We next 
note that T nm = T° m + 5T nm + O^ 1 ] where the Drude 
contribution, T^ m , is 0[N] and the weak localization con- 
tribution, ST nm , is 0[1]. If we now expand g for large N 
and collect all 0[l]-terms we get 



LR 



(T L D 3 ) 2 «R3 + (T R D 3 ) 2 ST L3 



PS 



R3J 



(7) 



For a cavity perfectly connected to all three leads (with 
Wl,r,3 S> Ap), the Drude and weak localization results 
for U = (at finite-Tg ) [l5| can be substituted into 
Eq. |(7J), immediately giving Eq. @ with f given by r^ 1 . 

To connect with the numerics of Ref. [X 31 ] . we now 
consider a tunnel-barrier with finite transparency < 
ps < 1 between the cavity and the dephasing lead. In- 
troducing tunnel-barriers into the trajectory-based the- 
ory of weak localization is detailed in Ref. [2l|. It re- 
quires three main changes to the theory in Ref. [l5j ]. 
(i) The dwell time (single path survival time) becomes 
t di = ( T o L )~ 1 I2 m PmW m . (ii) The paired path sur- 
vival time (for two paths closer than the lead width) 
is no longer equal to the dwell time, instead it is 
T D2 = ( T o^) _1 EmPmt 2 ~ Pm)W m because survival re- 
quires that both paths hitting a tunnel-barrier are re- 
flected (iii) The coherent-backscattering peak con- 
tributes to transmission as well as reflection, sec Fig. [3] 

For the Drude conductance we need only (i) above, giv- 
ing us T° TO = p m p n N m N n /K, where K = J2 k Pk N k- For 
the conventional weak localization correction we need (i) 
and (ii). The contribution's classical path stays within 
W of itself for a time T\y(e)/2 on either side of the en- 
counter (dashed region in Figs. O and [3]), thus we must 
use the paired-paths survival time, td2, for these parts 
of the path. Elsewhere the survival time is given by tdi- 
We follow the derivation in Ref. [HI with these new ingre- 
dients, and the conventional weak localization correction 
becomes 6T^l = -(p m p n N m N n /7<S 2 )(TDx/TD2) exp[-6]. 



The exponential with 9 = T£ P /td 2 + (rp; — r p P )/ T Di, is 
the probability that the path segments survive a time r^ p 
as a pair (t^ p Jl cither side of the encounter) and survive 
an additional time (r^ 1 — Tg P ) unpaired (to form a loop of 
length Tg 1 ). However, we must include point (iii) above 
and consider the jailed coherent-backscattering, shown in 
Fig. [3] We perform the backscattering calculation follow- 
ing Ref. [l5( (see also [l8j|) but using td2 when the paths 
are closer than W and Tdi elsewhere. We then multi- 
ply the result by the probability that the path reflects off 
lead to and then escapes through lead n. This gives a con- 
tribution, 5T^ = -{p m {l - p m ) Pn N m N n /N 2 ) exp[-6] 
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Figure 3: (color online) A failed coherent-backscattering con- 
tribution to conductance, 5T^. It involves paths which re- 
turn to close but anti-parallel to themselves at lead m, but 
reflects off the tunnel-barrier, remaining in the cavity to fi- 
nally escape via lead n. The cross-hatched region is when the 
two solid paths are paired (within W of each other). 



assuming n ^ m. There is a second such contribution 
with m <-> n. Summing the contributions for m =/= n 

ST nm = {p mPn N m N n /N 2 ){p m + p n -N/N)c- e (8) 

where N = J^kPk^k- If only the dcphasing lead has 
a tunnel-barrier, substituting the Drude and weak lo- 
calization results into Eq. |(7J), we obtain Eq. © with 
f = (1 — p)r^ p + t e ' . In this case the exponential in 
Eq. ([2|) is the probability that a path does not escape 
into the dcphasing lead in either the paired-region or the 
extra time (r^ 1 — r^ p ) unpaired (for the loop to form). 

To generalize our results to j dephasing leads, we 
expand the relevant conductance formula in powers of 
N and collect the 0[l]-terms. Then g wl = ST hR + 
Ei=i(^m | '^"i+^'n^Km), where the sum is over all de- 
phasing leads. The prefactors A m ,B m are combinations 
of Drude conductances and thus independent of Tg,Tg P , 
we need them to get power-law dephasing. However we 
can already see that there must be exponential decay 
with [(1 — p)t^ + Tg 1 J/t^, as for j = 1, where now t^ 1 = 
(roLy^mPmWm and pr- 1 = (toL)- 1 £ m p 2 m W m . 

Conclusions. We first observe that the dephasing- 
lead model has no independent parameter £. To our sur- 
prise it is the Fermi wavelength, not the dephasing-lead's 
width, which plays a role similar to £. Thus a dephasing- 
lead model cannot mimic a system-environment model 
with £ ~ L at finite . Our second observation is that 
Eq. @ with t = T£ is for a regime where U does not af- 
fect the momentum/energy of classical paths. Therefore, 
it does not contradict the result with f = r^ 1 in Ref. [l2| , 
valid for £ so small that dephasing occurs via a "single 
inelastic process with large energy transfer" [22[ ■ Intrigu- 
ingly their result is similar to ours for the dephasing-lead 
model. Could this be due to the destruction of classical 
determinism by the dephasing process in both cases? 

We finally note that conductance fluctuations (CFs) 
in the dephasing-lead model exhibit an exponential de- 
pendence oc exp[— 2te/t0] for p 3 <C 1 (HI, but recover 
the universal behavior of Eq. (fTJ) for p% = 1 (23|. How- 
ever external noise can lead to dcphasing of CFs with 



a TE-indcpendcnt exponential term [23|, similar to the 
one found above for weak-localization in the system- 
environment model. Thus our conclusion, that dephasing 
is system-dependent in the deep semiclassical limit, also 
applies to conductance fluctuations. 
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